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Abstract 

The Friedel oscillations in the vicinity of a Friedel-Anderson (FA) impurity are 
investigated numerically. For an FA impurity in the local moment limit the normalized 
amplitude A{(^) is S-shaped, approximately zero at short distances, approaching two at 
large distances and crossing the value one at the characteristic length ^1/2- Surprisingly, 
the Friedel oscillations of a simple non-interacting Friedel impurity with a narrow 
resonance at the Fermi level show a very similar behavior of their amplitude A(^). A 
comparison correlates the resonance width and the Kondo energy of the FA impurity 
with the characteristic length (^1^2 of the Friedel oscillations. 

PACS: 75.20.Hr, 71.23.An, 71.27.+a , 05.30.-d 

1 Introduction 

The properties of magnetic impurities in a metallic host were first studied theoretically by 
Friedel [Jj and Anderson |2j. Kondo [3J showed that a magnetic impurity with spin- flip 
scattering develops a singular behavior at low temperatures. As a consequence the ground 
state is a singlet state with zero effective magnetic moment. Schrieffer and Wolff [4j showed 
that the Friedel-Anderson (FA) Hamiltonian can be transformed into a Kondo Hamiltonian 
plus a number of additional terms. Its ground state is also a singlet state. The disappearance 
of the magnetic moment at low temperatures, the Kondo effect, is one of the most intensively 
studied problems in solid state physics , 0, [6], [7], [8], [9], [lO], [11], [12], [13], [11], [15], ^ 
[T7] . [T8] . [T9] . In the last decade the Kondo effect has experienced a renaissance. There 
is a growing interest in this field [20], extending from magnetic atoms on the surface of 
corrals [21j to carbon nanotubes [22], quantum dots [23], [21], [25], [26], [27], [28], [29] and 
nanostructures [30j. There are still many open questions, particularly the real-space form of 
the wave function and the resulting charge density and polarization. 



1 



Affleck, Borda and Saleur [31] (ABS) investigated the formation of Friedel oscillations in 
the vicinity of a Kondo impurity. The oscillating part of the electron charge has the following 
form. 



A (r) cos ( 2kFr + 



(1) 



where D is the dimension of the system, the coefficients Cn have the values Ci = 1/ (27r), 
C2 = 1/ (27r^) and C3 = 1/ (47r^) . They succeeded in calculating the asymptotic form of the 
function A (r) (which is a universal funtion of r/r^, tk = hvp/ {kBTK)= Kondo length) from 
first principles and performed a numerical calculation using NRG (numerical renormalization 
group). A{r) approaches the values two for r >> and zero for r << r^. One of the 
authors [32] applied the approximate solution of the FAIR theory (Friedel Artificially Inserted 
Resonance) to reproduce the result by ABS. The two calculations confirm each other. 

In the same paper one of the authors showed that a simple non-interacting Friedel im- 
purity with its resonance at the Fermi energy showed Friedel oscillations similar to the 
Kondo impurity. This reveals the interesting fact that the complex interacting Kondo impu- 
rity yields Friedel oscillations very similar to those of a simple resonance of a non-interacting 
electron gas. In both cases the amplitude A (r) approaches zero for r << tr and for r >> 
each spin contributes the value one to the amplitude. The transition happens at a charac- 
teristic length which is of the order of = Hvp/^ with vp being the Fermi velocity and F 
the half-width of the resonance. 

It had been observed earlier that the Friedel phase shifts and sum rules that were de- 
rived for non-interacting electron systems also apply to interacting electron systems (see for 
example [?]). This can be a great help in describing interacting electron systems, which are 
generally not very transparent. A very good example is the Kondo impurity with its Kondo 
resonance [?], [?], [?], [21J. However, the opinion of experts (in private conversations) are 
extremely divided, ranging from "there is really no Kondo resonance" to " yes there is a 
Kondo resonance but its form is very complicated and nobody can write it down". There- 
fore it is a considerable help when the complex system follows the same simple rules as a 
single electron system. In this paper we calculate the Friedel oscillations of the Friedel- 
Anderson (FA) impurity and demonstrate their similarity with the Friedel oscillations of a 
simple non-interacting Friedel impurity. 



2 Theoretical Background 

The FA Hamiltonian consists of spin- up and down free electron states c^^, a d- resonance dl 
with the energy Ed and s-d-hopping matrix elements VJ"^ between the conduction electrons 
and the d-impurity. A simultaneous occupation n^t and Udi of the d-resonance contributes 
a Coulomb exchange energy Uud-^Udi- 

{N-l N-1 ^ 

e^c^c^a + Vf[dlc^„ + cij^] + Eddid, \ + Uud^ndi (2) 
u=0 u=0 J 
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The FA impurity shows also the Kondo effect as Schrieffer and Wolff [1] have shown. The 
FA impurity has many degrees of freedom. Its behavior is determined by the d-resonance 
energy E^, the Coulomb energy U and the s-d-scattering matrix element IV^^^I . For simplicity 
we will discuss here the symmetric FA impurity with E^i = —U/2. 

For sufficiently large values of U/ \ Vsdf the FA impurity possesses a magnetic moment 
above the Kondo temperature. One difficulty in determining the magnetic moment is the 
fact that the magnetic state is not the ground state of the system. The ground state is a 
singlet state (or Kondo state) with vanishing the total moment, whose energy is lowered by 
the Kondo energy. For many realistic systems the Kondo energy is quite small. If warmed 
above the Kondo temperature the FA impurity is magnetic, but it is still at rather low 
temperature. We denote this as the magnetic pseudo-ground state. In the FAIR approach 
the magnetic wave function has the form 

n — 1 n— 1 



Aal^bl + Bal^d\ + Cd\hl^ + Dd\d\\ \{ aj^ J] ^4*^0 (3) 

i=l i=l 

Here aj^ and are two artificially inserted Friedel resonances with a composition 

V 

(As discussed in previous papers [33] the FAIR states have a relatively simple interpretation 
in Hilbert space.) The remaining states in the spin up basis are jctj-f^l which are orthogonal 
to each other and to aj^.ln addition the free electron Hamiltonian i^ot foi' spin- up electrons is 

\t |7T_.|^t 



made sub-diagonal in this basis, i.e. all matrix elements <^a]^^ l-f^otl '^/^t/ '^'^.a* z^, /i 7^ 0. As 
a consequence the whole basis is fully determined by aj,^. The FAIR state fojj^ with its basis 
1 1 has analogous properties. The details of the FAIR theory are described in previous 

papers [51], [33], |2S], etc. 

The mean field theory produces two d-resonances in the magnetic state, one for the spin- 
up band at Eii + Un^i and another one for the spin-down band at Eii + Uridf. Here Ud^, 
are the partial occupations of the spin-up and -down d-resonances. In mean field theory 
these resonances have the half width F = tt |\4rf|^p (p=density of states). Mean field theory 
completely decouples the s-d-scattering between different spin bands. In reality a doubly 
occupied d-state decays in both spin bands, and the decay rate is twice as large [37]. The 
FAIR theory of the magnetic state yields d-resonances with the correct resonance width [38J. 

The ground state of the FA impurity is a singlet state. The latter is obtained in FAIR 
by reversing all spins in equ. (E]). After ordering the spin sequences the two states are 
added and normalized. For the singlet state the composition of the FAIR states aj, and 
60 has to be optimized again (and is very different from the composition in the magnetic 
state). The two (opposite) magnetic states in @ are not orthogonal to each other and yield 
a finite interference between the two magnetic components. The strength of this interference 
determines the Kondo energy. 
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Dd\dl 



n: 



=1 "«t 



n 



A'6S,aS, + 5'4aS, + C%dl + D'd\dl U^^i K^i 



In the numerical calculation we use the Wilson band [12] which is an electron-hole sym- 
metric band, normalized by the Fermi energy and extending from (—1 < e < +1) with a 
constant density of states p = 1/2. The electron states in the Wilson band are represented 
by a finite number of Wilson states. We generally start our numerical calculation with 
= 60 Wilson states with energies of±|,±|,±^,.. ±230-, ±23^- The ratio between neigh- 
boring energies is A = 2 (except for the two states next to the Fermi level). The wave 
function of the Wilson states in real space is described in previous papers |39], [32]. From 

the composition of the two bases ^^4^ terms of the Wilson states cj, one obtains 

their wave functions in real space. 

The calculation in real space can be done in one, two or three dimensions. The results 
differ by the dimensional factor C^i/r^and the phase shift — -D| (see equ.([I])). If one splits off 
the dimensional factor then we expect that the amplitude A (r) is the same for all dimensions. 
Therefore we choose the simplest case, the one-dimensional FA impurity. 

The calculation of the real space density at the position r is in principle straight forward 
but rather tedious due to the interference terms between the different bases. Here one 
has to calculate a large number of multi-electron scalar products in form of determinants. 
Fortunately this has to be done only once and can be used for all distances. 



3 Numerical Results and Discussion 

We performed the calculation of the Friedel oscillation for different combinations of parame- 
ters. In a first series we used for the s-d-hopping matrix element the value IV^^^I = 0.03 and 
for the Coulomb energy the values f/ = 0.1, 0.2, 0.4, 0.6, 0.8 and 1.0 with Ed = -U/2. In 
the first round 60 Wilson states with energies of±|,±|,±^,.. ±2io,23o represented the free 
electron basis. As ABS already noticed the Wilson states with an energy ratio of A = 2 are so 
far apart in their energy that the results are not yet very accurate. (For example the Wilson 
state with the energy | is composed of all states in the energy range of i < e < l.This means 
that it has roughly an energy uncertainty of |. In simple words, the energy uncertainty of 
the Wilson states divided by the energy AEw/Ew is about 1/3 for A = 2). In the FAIR 
approach it is relatively easy to sub-divide the Wilson states geometrically [32j. The FAIR 
states aj and 6q can be well interpolated during this procedure. From the FAIR states one 

obtains the whole bases jc^lj and 

After the first sub-division we have a new Wilson basis of 120 states with an energy ratio 
of A = 2^/2_ Ti^g second sub-division yields 240 states with A = 2^/^ ^ 1. 19. Now even the 
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large energy levels lie close together, and the energy uncertainty /S.Ew I Ew takes the value 
of about 0.09. 

In the calculation we measure the distance from the impurity in units of half the Fermi 
wave length \p, i.e. ^ = 2r/\p. In these units the Friedel oscillations have the period 
"one". We multiply the actually calculated amplitude by the factor 2ttS, to cancel the one- 
dimensional prefactor in equ. ([1]). The resulting normalized amplitude v4(^) is then valid 
for all dimensions of the sample. 

In Fig.l the Friedel oscillation amplitude A{E,) is shown as the function of the logarithm 
of the distance £ = logg (0 for the different parameters. 
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Fig.l: The amplitude A {€} of the Friedel oscillations for a symmetric 
FA impurity with l^^^l = 0.03 and different values of U with = —U/2. 
The abscissa is the distance £ = log2 (0 from the impurity where ^ = 2r/\F- 

The right curve (stars) in Fig.l has the largest U = 1 and behaves very similar to a 
Kondo impurity that is discussed in ref. [31], [32]. For large distances the amplitude A{^) 
reaches the value 2, and for short distances the amplitude essentially vanishes. The transition 
between these two regions occurs at the £i/2-point where A {£) = 1. For U = 0.8 one obtains 
a similar curve as for U = 1.0, only shifted to the left. With decreasing U the £i/2-point for 
the amplitude moves to shorter distances 

For U = 1 and U = 0.8 the positions of the £i/2-points are at 9.4 and 13.1. If one shifts 
one of the curves by 3. 7 then one obtains a single curve. It represents the universal curve 
for the FA impurity with a full magnetic moment. If U is further decreased then the curves 
for A (^) are further shifted to the left and for U < 0.4 they don't reach the horizontal axis 
anymore. Essentially the whole environment of the impurity is filled with Friedel oscillations. 
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Surprisingly the curves develop a small minimum at about £ = 20 with decreasing value of 
U . 

For the smallest value of f/ = 0.1 one is in the perturbative region while for U = 1.0 
one is in the local moment region (see for example Krishna-murtha et al. [ID]). Generally 
the properties of the symmetric FA impurity are characterized by the parameter U/ (ttF) = 
U I (vr^ p) ■ In table I the values of f// (ttF) are collected for the different examples. Our 
curves for the Friedel oscillations approach the local moment behavior somewhere above 
Uj (ttF) ^ 4. 



u 




U/ttF 


r 


S 


0.1 


0.03 


0.66 


0.51 


0.303 


0.2 


0.03 


1.33 


0.67 


0.354 


0.4 


0.03 


2.67 


0.83 


0.434 


0.6 


0.03 


4.05 


0.90 


0.479 


0.8 


0.03 


5.33 


0.93 


0.497 


1.0 


0.03 


6.67 


0.95 


0.500 



Table I: The effective coupling strength f/ / (ttF) for the different 
parameter of Fig.l. r and 5* are discussed in the text. 

3.1 The magnetic half of the singlet state 

The FAIR solution for the FA-impurity consists of two entangled magnetic states as shown 
in equ. ([1]). Each half in equ. 01]) possesses a magnetic moment. For large U j (ttF) the 
magnetic moment is relatively well defined because the sum S = {A^ + 5^ + + D'^) is 
close to 0.5 and almost normalized to 1/2. If we define r = {B'^ - C^) / {A^ + B"^ + C'^ + D'^) 
then the moment has the value ix/ixb ~ r as long as ~ 0.5. However, with decreasing 
U/ (ttF) the sum S becomes considerably smaller than 1/2 because it is the singlet state 
that is normalized, and the interference between the two magnetic parts becomes more and 
more important. The values of r and S are also collected in table I. It should, however, be 
emphasized that the magnetic moment of the magnetic half in the singlet state is not the 
same as the magnetic moment of the corresponding magnetic state as presented by equ. ([3]). 
Although the structure of the states is the same, the FAIR states a\ and h\ assume rather 
different composition. 

For U = l,Ed = -0.5 and IKdl = 0-03 the magnetic halves of the singlet state are very 
well normalized to the value of 5 = 0.500. Therefore it is interesting to study their properties. 
In Fig. 2 the half of the ground state with the net d-spin down is investigated. Both the spin- 
up and the spin-down bands of the conduction electron gas show charge oscillations with 
the period one. Their magnitudes are identical within the accuracy of the calculation. At 
large distances they reach the value A = 0.5. This means that the total charge oscillation of 
both spins for both magnetic components yields the value two which agrees with the results 
of Fig.l. However, for short distances one observes for each spin band a value of 0.05. This 
value appears to contradict the vanishing Friedel oscillations in Fig.l. 



6 



0.5- 

< 0.4- 

"D 

1. 0.3 ^ 

E 

CO 

>^ 

■■^ 0.2- 

c 

<a 

T3 

0.1- 
0.0- 



magnetic part 

N=4*60 
geometric 



/ 



— ^— spin up 
—V— spin down 



IV, 



0.03 
Eh = -0.5 

d 

u = 1 



FO R3V03En50U100 b 



— I— 
5 



— I — 
10 



15 20 

^ = iog,(^) 



— I 
25 



Fig. 2: The amplitude A{^) of the charge oscillations of the spin- up (full 
triangle up) and the spin-down (open triangle down) electrons in the 
magnetic component of the singlet state. 

The reason for this apparent contradiction hes in the phase of the different oscillations. 
The position of the maxima of the charge oscillations may be at ^ma- In the absence of a 
phase shift the maxima of the Friedcl oscillations would be at integer values of ^. Therefore 
the deviation of the from an integer, S^, yields the phase shift in units of 2tt. The 
positions of 5^of the maxima and minima of the charge oscillations of the spin-up part and 
the spin-down part of the conduction electrons are plotted in Fig.3. The full triangles show 
5$,mi of the minima of the charge oscillations and the open triangles represent the maxima. 
In addition the up-triangles represent the spin-up electrons and the down-triangles the spin- 
down electrons. 

For large distances the maxima positions for both spin orientations he at 0.75 and the 
minima at 0.25. Therefore the charge oscillations add constructively in this region. However 
for distances which are smaller than ^1/2 the maxima and minima of spin- up and down 
electrons are shifted by ±0.25 yielding a relative difference of 6C, = 0.5, corresponding to a 
phase difference of tt. Therefore the charge oscillations of spin-up and down electrons cancel 
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in this region. 
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Fig.3: The phase relations of the charge oscillations of the magnetic half of 
the FA impurity with U = 1, Ed = —0.5 and \Vsd\ = 0.03 as a function of the 
distance i = log2 {^). The full and open triangles show the positions (modulo one) of the 
minima and maxima of the Friedel oscillations for the spin-up and spin-down 
electron (up and down triangles). 

While the total charge oscillation in the magnetic half of the singlet solution cancels 
essentially to zero in the region of i << ii/2 one obtains a polarization oscillation in this 
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region. In Fig.4 the amplitude of this polarization in plotted as a function of £ — log2 {^). 
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Fig.4: The amplitude of the polarization oscillation of the 
magnetic magnetic component for the U — 1 FA impurity 
as a function of the distance £ — log2 

If one performs the same analysis for the FA impurity with U = 0.1, = —0.05 and 
\Vsd\ = 0.03 one observes essentially two important differences. Fig. 5 shows that (i) the 
phase difference between spin-up and down electrons at short distances is much less developed 
than in the local moment example in Fig. 3, and (ii) the transitional region where the phases 
split is shifted to much shorter distances. The fact that the phase differences at short 
distances is much less than tt is the reason that the two charge oscillations of spin-up and 
down electrons do not cancel. The Friedel oscillations maintain a relatively large amplitude 
even at short distances. 
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Fig.5: The phase relations of the charge oscillations of the magnetic half of 
the FA impurity with U = 0.1, = —0.05 and = 0.03 as a function of the 

distance ^. The full and open triangles show the positions (modulo one) of the 
minima and maxima of the Friedel oscillations for the spin-up and spin-down 
electron (up and down triangles). 

3.2 Fermi liquid picture 

Nozieres [1T| derived from Wilson's renormalization fT2] and Anderson's scaling theory |12] a 
Fermi liquid description of the Kondo effect at low temperatures. In this description, which 
should also apply to the FA impurity in the local moment regime, the magnetic moment 
forms a strongly bound singlet state with a conduction electron. This removes the magnetic 
moment as well as one conduction electron from the system. Virtual excitations cause a 
weak interacton between the quasi-particles. In ref |32] one of the authors pointed out the 
similarity between the Friedel oscillations of a Kondo impurity and a Friedel impurity with 
a very narrow resonance at the Fermi level. This is another confirmation of the Fermi liquid 
description of a Kondo impurity. In the present paper we want to show that this simililarity 
also extends to the FA impurity. 

Fig. 6 shows the Friedel oscillations of a simple Friedel impurity with spin- up and down 
sub-bands. The resonance d-state lies at the Fermi level, i.e. = 0. The s-d-hopping 
matrix element \Vsdf is varied by four orders of magnitude between 10^^ and 10^^. 

There is a remarkable similarity between the Friedel oscillations of the FA impurity in 
Fig.l and of the Friedel impurity in Fig. 6. If we consider first the Friedel impurity with 
\Vsdf = 10~^ then the overall shape of the A {£) curve is very similar to the A (£)-curve for 
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the FA impurity with U — 1. Both curves approach the value two at large distances and 
zero at short distances. For the Friedel impurity an increase in \Vsdf by a factor of 10 shifts 
the curve by log2 (10) to the left. The two left curves in Fig.6 no longer reach the value zero. 
This must be due to the fact that F is no longer sufficiently small compared to the band 
width of one. 

We observed a similar behavior for the FA impurity when we reduced the value of U. The 
similarity with the Friedel impurity suggests that the reason here is also due to the finite 
band width and not due to the fact that the FA impurity with U / (ttF) < 3 is no longer in 
the local moment limit. 

For the Friedel impurity with \Vsd\^ = 10^^ the position of the £1/2 point is at £1/2 ~ 12.6, 
corresponding to a length of ^1/2 = 2^^-^ fa 6208. The product between ^1/2 and the resonance 
half-width F = tt \Vsdf | = 1. 57 x 10"^ is ^1.57x lO'^ * 2^2.6 g.lO (F is measured 
in units of the Fermi energy Ep and ^ in units of half the Fermi wave length \p/2). This 
product is universal (as long as F is sufficiently small compared with band width of the 
system) . 
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Fig.6: The amphtudes of the Friedel oscillations for different 

I 1 2 

Friedel resonances with Ea — O and \Vsd\ = lO"-' with j — 1,2, 3, 4, 5. 

We want to quantify this similarity. For this purpose we determine for each curve in Fig.l 
and Fig.6 the distance ii/2 — log2 {^1/2) where the amplitude A {£) has half the saturation 
value A {£1/2) — 1- For the Friedel impurity we plot log2 F of the resonance half- width F as 
a function of £1/2 = log2^i/2- That yields the full points in Fig. 7. 

The similarity between the Friedel curves in Fig.6 and the FA curves in Fig.l supports 
the interpretation that the FA impurity has a resonance at the Fermi energy. However, its 



11 



half-width is not well known. We assume that it is proportional to the Kondo energy (which 
is differently defined in different theoretical approaches). We choose for the FA impurity 
in Fig.l the definition of the Kondo energy by means of the susceptibility. The latter we 
obtain through a linear response calculation [13] yielding Ey.. In Fig. 7 we plot the logarithm 
of this Kondo energy log2 (E^) versus £1/2 = log2 {^1/2} for each curve in Fig.l (stars). Both 
sets of points for the FA and the Friedel impurity (circles and stars) lie on two straight lines 
with the slope of —1 and are separated by a vertical distance of 1.25 which corresponds to a 
factor of 2.4. If we assign to the FA impurity a resonance half-width of TpA ~ ^x/'^-^ then 
the values for Tp and F^^ fall on the same straight line with the slope of —1. 

Our calculations of the Friedel oscillations not only show that there is a Kondo resonance 
at the Fermi level but also suggest that even for relatively small U / (vrF) the quasi-particle 
behavior is similar to that of a Friedel impurity. 




5 10 15 



Fig. 7: The logarithm of the half- width Tp for the Friedel impurity (full circles) 
and the logarithm of the Kondo energy for the FA impurity are plotted 
versus the corresponding values £1/2 = log2 (^1/2) -values. The latter are obtained 
from Fig. 6 and Fig.l as the position where A (£) takes the value A (^1/2) = 1- The 
two full lines are straight lines with the slope minus one. Their vertical 
separation is 1.25, corresponding to a ratio of 2^'^^ ps; 2. 4 between Kondo 
energy and Kondo resonance half width. 

4 Conclusion 

In this paper we investigate the Friedel oscillations of a set of symmetric Friedel- Anderson 
impurities. In the local moment limit (large U/irT) the normalized amplitude A{^) essen- 
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tially vanishes for short distances and assumes the value 2 for large distances. In this range 
oi U / (ttF) the A (^)-curves show universal behavior because curves for different U / (ttF) 
can be shifted into a perfect overlap. The physics behind this behavior of the amplitude 
of the Friedel oscillation is illuminated by the study of the simple Friedel impurity with a 
very narrow resonance at the Fermi level. For these non-interacting Friedel impurities one 
obtains real space Friedel oscillations which are very similar to those of the interacting FA 
impurities. 

In the case of the Friedel impurity one can derive two conclusions from the Friedel 
oscillations: 

• The amplitude A{^) reaches the saturation value of two only when the resonance lies 
exactly at the Fermi level. 

• The .^i/2-point where the normalized amplitude A (C1/2) is equal to one yields the 
resonance half- width Tp through the universal relation VFii/2 ~ 0.10. 

The similarity between the Friedel oscillations of the FA and the Friedel impurity supports 
the concept of a " Kondo" resonance at the Fermi level. The position of the ^1/2-point suggests 
that the half- width of the FA resonance is of the order of the Kondo energy F^^ ^ Ey,/2A . 

Our calculations of the Friedel oscillations support the concept that there is a Kondo 
resonance at the Fermi level with the phase shift 7r/2 and suggests that even for relatively 
small values of U / (ttF) the quasi-particle behavior is similar to that of a Friedel impurity. 

It should be emphasized that the spatial behavior of the FAIR wave function and its 
charge and spin oscillations provides quantitative information about Kondo energy, Kondo 
length and resonance width of a FA impurity without the need of a magnetic field or excited 
states. 
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